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Abstract 
Introduction. Domestic and foreign works on the milling of complex-profile surfaces with a ball-end tool were 


analyzed. Methods of surface quality control and ways to provide amplitude parameters of roughness, based on research 
data and field experiments, were considered. Theoretical provisions on the determination of cutting forces and the 
results of vibroacoustic diagnostics were presented. 

Materials and Methods. The methods of correlation analysis, comparison and generalization of the results were 
applied. The data were calculated at different tool angles, taking into account the instantaneous cutting forces, and were 
fixed in the range of values of the variable feed per tooth (fz) and the angle of inclination of the surface (y). The 
vibroacoustic diagnostic data and theoretical data of the presented model at different tool inclination angles were 
verified by experiment. Consequently, such methods can be used to predict surface roughness parameters. 

Results. The relationship between cutting forces, tool inclination angle, and vibroacoustic diagnostics data was found. A 
model of the cutting force and tool displacements was formulated taking into account the inclination of the surface. The 
optimal range of the inclination angle of the tool to the surface to be machined, at which the minimum values of the 
amplitude parameters of roughness were achieved, was determined. The sound vibrations obtained empirically, 
presented in spectral and wave forms, were in good agreement with data from other sources. This allowed us to 
conclude about the feasibility of forecasting and monitoring roughness parameters in real time through acoustics. 
Discussion and Conclusions. It was established that the growth of forces in the direction a,(X) and fz(Y) was observed 
at y > 40°. This was due to the distribution of the components of the cutting force along the cutting edge and depended 
on the inclination of the surface. The amplitude parameters decreased when the angle increased from 10 to 40 degrees. 
The found interrelations of force analysis, processing directions, and vibroacoustic diagnostics have validated the use of 
vibroacoustic diagnostics to predict surface roughness. Acoustic diagnostics, regardless of the layout of technological 


equipment, enables to quickly adjust the sound device and assess the impact of cutting modes on roughness parameters. 
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Introduction. Specialists do not have enough software tools that could help predicting the value of the cutting force, 
vibration amplitude, etc., when programming machines with numerical control (CNC). As a result, it is required to 
select the feed rate based on experience, setting a constant feed rate for the entire processed complex surface. This can 
reduce the processing efficiency and surface quality. To solve this problem, it is required to build a highly accurate and 
reliable model for predicting surface roughness. 

Currently, two main approaches are used in such forecasting. In the first case, photoelectric sensors transmit 
roughness data online. This method efficiency is quite low because chips are formed during the cutting process, and a 
lubricating fluid (lubricoolant) is supplied to the area of the recorded microrelief. The second group includes control 
methods using artificial intelligence, as well as theoretical modeling and modeling of empirical regression analysis. In 
practice, methods from the second group are most common. 

The authors of paper [1] explored the spherical tool milling of aerospace components made of aluminum alloy LM6 
and developed a model using ARMAX structures (autoregressive moving average with exogenous inputs) to predict 
surface roughness under various processing modes. In [2], a mathematical model of the trajectory of a spherocylindrical 
tool for processing mushroom-shaped blades of a steam turbine was presented to improve the processing efficiency and 
surface quality. The authors [3] created an analytical model of the milling with a spherocylindrical tool and considered 
the situation of vibration occurrence to predict the limits of tool stability at different angles of inclination and spindle 
rotation speed. 

In [4], the processing of the impeller on a five-axis CNC machine using a spherical cylindrical tool was studied. The 
tool trajectory and optimal cutting modes had to be developed due to the complex curved surface and the close location 
of the wheel blades, since in this case, it was impossible to use standard CAM systems (computer-aided manufacturing - 
an automated system or its module for the preparation of control production programs). 

In [5], a system for optimizing the feed rate for three-axis milling with a spherical tool was described. The authors 
used the geometric modeling capabilities of B-rep ACIS: miter cutting was converted to an orthogonal model for 
various tools and workpieces. The authors [6] created a mathematical model of a constant cutting force [7] for 
processing a complex curved surface with a spherocylindrical tool on a CNC machine through maintaining a constant 
cutting speed at different tool inclination angles. Paper [8] presents a dynamic model of the milling with a spherical 
cylindrical tool, including cutting modes and tool edge wear. The authors [9] developed a real-time monitoring system 
with error compensation to improve accuracy in the production of free-form parts. In [10], the energy consumption of 
various cutting strategies was studied to evaluate the efficiency of each cutting strategy. The authors of this work have 
established that the optimal strategy enables to provide the specified parameters of surface quality. Many scientists used 
the Taguchi method to optimize processing parameters [11—13]. Some of the surveys are narrowly focused, and more 
work is needed to scale up the results [14-18]. 

The analysis of literature sources found the absence of recommendations on using the method of forecasting and 
monitoring roughness in real time. This is even more so in cases where sound vibrations in the air are taken into 
account. 


Materials and Methods. To determine the cutting forces, a model from [7, 19] is used: 


dF =Kydl, + KAS, 
dF, =Kyedl; +KyedSz, 


(1) 
dF, =Kgedl; + KgcAS.;. 


Here, Kj, Kye, Kae — edge specific coefficients; K,,, K;-, Kae — shear specific coefficients; dl; — incremental cutting- 
edge length /; [20]; S.; — cross-sectional area of the cut segment. 

In this model, a set of curvilinear coordinate systems, perpendicular to the tangents to the spherical surface, is used 
to set the resulting force acting on the i-th infinitesimal segment of the cutting edge. 


The geometry of the spherocylindrical end mill is schematically shown in Figure 1. 
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Fig. 1. Diagram for determining: a — cutting forces of a spherocylindrical cutter; 


b — segment to be cut; c — cross-sectional area of the given figure [16] 


The diagram shows the angle of inclination 4 of a milling cutter with a spherocylindrical initial surface of the tool 
along the cutting edge. Here, we also see: 

— tangential dF, radial dF, axial dF; cutting forces acting on the j-th tooth, and the coordinates of the tool for the 
spherocylindrical end mill (Fig. 1 a); 

— chord L,,;,, segment angle ®,, arc length /,,, (Fig. 1 b); 

— cross-section area (Fig. | b, c). 
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= . ht 2 
Sj = Lyon * Maver: @) 


where L>,;, — length of the second chord, mm; h’,,,., — average chip thickness, mm. 


(3) 
I very = «SZ + SIMA. (4) 
Length of the first chord L;,;, arc length /,,,, and circle segment angle ®;: 
L, = 2-sin- lwe . R, (5) 
ch 2R 
®O,-R 
lore = > 
ze (6) 
@, = 2arccos [1 -2- =). 
2R (7) 
Angle A (Fig. | c) for calculating the average chip thickness: 
A = arctan [4] : 
Lin / 2 (8) 
Radius vector r of the current point of the initial tool surface: 
Rsing cosQ. 
Rsing sinQ 
ra aly (9) 
Reos@ 
1 


The equation of the cutting edge is solved as the equation of a curve located on the original instrumental surface and 
crossing the meridians at angle 190°. A parametric equation of the form © = () describes a certain line on the 
sphere. Cosines of the angles between the axes of XYZ coordinate system of the spherocylindrical cutter and the tangent 
line to the curve Q = Q (q) on the initial tool surface: 

dX cos@ cosQ do — sing sinQ do 


cose =, 
dS Jde + sin’ dQ? 


dY cos@ sinQ d@ — sing cosQ do 


cosB = a P 
dS Vdo + sin dQ? (10) 


dZ — sing d@ 


cos6 = = 7 
dS do + sin“@p dQ 


where dS — cutting edge arc differential. 
Cosines of the angles between the coordinate axes and the tangent to the meridian (at Q = const) on the initial tool 


surface are represented as Radzevich’s dependences [21]: 


dX 
cosa = —— = cos@ cosQ, 
dS 


cospB = — = cos@ sinQ, (11) 


dZ ; 
cos6 = — = -— sing. 


ds 
We set the condition: Q = i. Then, the cutting-edge equations in parametric form: 
Reosi 
chq (2% + C)- 
Rsinx (12) 
chq (A + C) 


Z = Rthg(r% + C), 
where ch — hyperbolic cosine, th — hyperbolic tangent. 
To determine the infinitesimal length of the cutting edge, we can use the expression proposed in [7, 16, 19]: 
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dl = sr’ (w) + (r' (wy) + (<' (wy) ay, (13) 
where 
r(y) = Vl - (y cot a — 1) 
toe —R (y cot A — 1) cotr . (14) 
V1 - (wy cot A — 1) 
z’(w) = X* + ¥7 + (R- ZY cota 
Main Text 


Experiment Details. At the preparatory stage, the literature was studied [22-27]. The data for planning the 
experiment are described below. Machining was carried out in passing and counter directions (Fig. 2), on blanks with 
the properties of steel 45 (foreign analogues — C45, 1045). We used a spherocylindrical end carbide milling cutter with 
a diameter of D = 8 mm, with two teeth from Sandvik Coromant (z = 2). The ratio of the overhang of the tool installed 
in the chuck was taken as | / D, equal to 4. The milling tool made of fine-grained tungsten carbide had anti-wear coating 
TiAIN and parameter 4 = 30°. Feed per tooth fz = 0.2 mm/tooth, allowance for all samples a, = H = 0.2 mm, lateral 
pitch a, = 0.2 mm. 

To provide equal cutting speed V, at different angles of inclination (y) of the tool to the surface to be processed, the 
rotation frequency (n) varied in the range from 1478 min ' to 8000 min‘. 


Z 
a, 
~ el 
a) 


Fig. 2. Types of the conditional contact surface of workpiece and tool at different values of tool inclination angles in the direction: 


b) 


a — passing; b — counter 


Angle y was selected from the range from 10° to 50°; therefore, the calculation of the effective diameter Deap (Lich) 
was carried out using the formulas: 


= 2Dcosm = 2. R - cos (90° — y) 
y > O° 40r (15) 
l 
D 7 = 2-sin:-—.-R. 
cap Lich OR 

In the experiments, a DMG DMU 50 CNC machine (Heidenhain control system) was used. The surface roughness 
after treatment was measured using profilometer-profilograph-contorograph Surfcam 1800D. In accordance with the 
current standards', measurements were carried out along various routes. To find the maximum values of the surface 
roughness parameters, 50 % Gaussian filtration was used. Vibroacoustic diagnostics was performed according to GOST 
R ISO 7919-3-997 using spectrum analyzer Zet 017-U2. As an output evaluation of the processing efficiency, amplitude 

A of the vibrations of the spindle assembly (mm), varying with time (tf, s) was used. 


' GOST 27964-88. Measurement of surface roughness parameters. Terms and definitions. Moscow, 2018. 13 p. (In Russ.) 
? GOST R ISO 7919-3-99. Mechanical vibration of non-reciprocating machines. Measurements on rotating shafts and evaluation criteria. Moscow, 


2000. 12 p. (In Russ.) 
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Research Results. The experiment has shown that the real surface roughness differs significantly from the data of 
theoretical models — the results of the kinematic-geometric projection of the tool on the workpiece. Currently, there is 
a growing demand for details of a spatially complex configuration and surface quality requirements. In this regard, it is 
worth focusing on the requirements of predictability of machining results, expressed, among other things, by roughness 
indicators. 

Tables 1 and 2 present the results of the experiment of amplitude parameters of roughness described in 
GOST R ISO 4288”. 


Table 1 
Roughness parameters depending on the machining direction 
with a spherocylindrical tool under climb milling 
Amplitude parameters of roughness, um 
Angle, ° Ra Rq Rz Rt Rp Rv 
10 0.457 0.569 2.912 5.000 1.671 1.240 
20 0.597 0.735 3.936 5.159 2.335 1.600 
30 0.449 0.561 3.007 4.960 1.808 1.200 
40 0.479 0.590 3.232 4.239 1.823 1.407 
50 0.787 1.027 6.615 12.68 4.119 2.496 
Table 2 
Roughness parameters depending on the machining direction 
with a spherocylindrical tool under up-cut milling 
Amplitude parameters of roughness, 1m 
Angle, ° Ra Rq Rz Rt Rp Rv 
10 0.806 1.006 4.776 5.960 2.776 2.000 
20 0.752 0.935 5.207 6.440 2.847 2.359 
30 0.509 0.647 3.727 5.239 2.223 1.503 
40 0.628 0.805 4.832 7.400 2.912 1.919 
50 0.730 0.915 5.111 7.639 3.039 2.072 


Analyzing the data of Tables | and 2, we note that the prevailing influence on the “height” of the roughness is 
exerted by the cutting direction. The quantitative value of the amplitude parameters is “lower” under climb milling. 
Hereafter, we will talk only about this type of milling. Based on the data in Table 1, a correlation analysis was 
performed, and a correlogram was made up (Fig. 3). The found dependences were correlated with the results of 
papers [28, 29]. 


— 0.9...1.00 0.8...0.89 --~  0.7...0.79 


Fig. 3. Correlogram of the amplitude roughness parameters after milling under climb milling 


The equations obtained during the correlation analysis helped us to move away from the direct normalization of 
roughness parameters. At the same time, analytical dependences provided minimizing the number of controlled factors. 
This reduced the complexity of the transition to real-time monitoring of sound vibrations. 


3GOST R ISO 4287-2014. Geometrical Product Specifications (GPS). Surface texture. Profile method. Terms, definitions and surface texture 
parameters. Moscow, 2015. 20 p. (In Russ.) 
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In the presented study, the cutting forces (Fx, Fy, Fz) were calculated (Fig. 4). It can be seen that the inclination 
angle of the surface y affected significantly the cutting forces. 


Fi,N Fi, N 
100 100 
—Fy Fy 
50 Fz 50 —F, 
Fx Fx 
0 0 
—50 50 
—~100 i 100 ts 
0 0.005 0.01 0.015 0.02 t, 0 0.005 0.01 0.015 0.02 
a) b) 


Fig. 4. Cutting forces calculated for: a— y = 10; b— y = 40 


The data obtained allowed us to establish that the forces decreased with an increase in the tool inclination angle 
from 10° to 40°. In addition, force Fx was most sensitive when the inclination angle of the surface changed. It should 
also be noted that the cutting forces calculated according to the above-mentioned method were in good agreement with 
the empirical results obtained [6, 17, 24, 26, 27, 30]. 

With an increase in angle y, the maximum area of plunging of the tool into the workpiece decreased (determined by 
the working angle and the active number of teeth z,). Consequently, when finishing milling with a spherocylindrical 
tool, “pulsating” forces may be present, since for the inclination angle of the surface y > 0°, the number of active teeth 
is often less than one (z, < 1). 

To establish z,, sound vibrations during cutting were analyzed. The authors of [31] proposed to divide the spectrum 
of vibrations of working bodies into ranges: 

— 20...300 Hz — low frequency, 

— 300...1500 Hz — medium frequency, 

— 1500 Hz and more — high frequency. 

The experimentally obtained range of sound vibrations is shown in Figure 5. 


Fig. 5. Experimentally obtained spectral type of sound wave (segment is | sec) 


The vibrations of the mounting elements and feed drives should be attributed to the low frequency range, the 
frequency of tool, spindle, etc. — to the medium and high frequency ranges. In this work, from the point of view of 
acoustic diagnostics, optimal cutting modes were determined that supported the stable operation of the technological 


system. 
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The frequency corresponding to the cutting frequency was pre-determined (Fig. 6). Filtering was applied to sound 
vibrations under: 

— spindle assembly movement (Y, Z), 

— rotation (n), 

— table movement (X), 


— noise in the laboratory. 


Fig. 6. Spectrogram of the sound wave after filtering (segment is equal to 1 sec) 


Some spread in values (Fig. 6), associated with a huge number of factors that were difficult to exclude when 
conducting a full-scale experiment, is noteworthy. 

The middle of the pronounced frequency range shown in Figure 5 corresponded to a rotation of 1500 rpm with a tool 
tilt of 40°. Acoustic diagnostics confirmed that the cutting went on almost continuously, i.e., without shock “pulsating” 


loads. The latter situation is shown in Figure 7. Here, the allocated evaluation time is 0.5 s, the number of periods is 25. 


Fig. 7. Filtered fragment of the sound wave based on the sound of cutting at y = 40° (segment is 1/2 sec) 


Figure 7 shows an ambiguously interpreted sound wave. However, it is quite easy to characterize it if we delve into 
the essence of mechanical milling. Under the climb milling and rounding of the cutting edge, plastic deformations of the 


metal with residual stresses remain on the treated surface. Even with proper coolant supply, temperature deformations 
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occur in the surface layers of the material and distort the crystal lattice (residual stresses of the HI kind). Along with 
this, after removing the tool, the stretched upper layers of the metal acquire residual compression stresses. If we 
describe the machining from the kinematics of the process, of note, the chip removal factor, namely value ©, — the 


cutting angle of the material with the inclined spindle in the passing direction (Fig. 8). 123 
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Fig. 8. Segment 2 to be cut off during machining: a — top view; b — projection of the segment being cut [16] 


We specify that the tool cut in zone 2. The data obtained were validated by the simulation results, but were 
considered for the ideal case. 

Along with these factors, theoretical and empirical discrepancies were often caused by vibrations associated with 
tool deflection and wear, as well as with the geometry of the cutting zone. Figure 7 shows the process of cutting the 
workpiece material with the j-tooth of the tool, which, along with other factors, determines such an interesting shape of 
the resulting wave. 

In addition to the above, Figure 9 presents dependences which show that with a decrease in the amplitude of 


vibrations, the angle of inclination increased from 10° to 40°. 
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Fig. 9. Vibration amplitudes depending on the surface inclination angle: solid line — 10°, dotted line — 40° 


M. R. Gimadeev and V. V. Gusliakov* analyzed the international practice of profile assessment of surface 
roughness. The authors of work [32] put forward proposals for the transition to a three-dimensional assessment. In the 
presented study, roughness measurements were carried out in an orthogonal coordinate system to make a 
comprehensive assessment. 

The analysis of the profilograms obtained (Fig. 10) allowed us to conclude about the regularity of the surface 
microrelief after processing with a spherocylindrical end mill. In this case, it is impractical to take into account the 
roughness of the tool, since its manifestations are insignificant. Some additional increase in roughness was observed in 
the direction of the side pitch a,, which verified the increase in elastic squeezes of the tool. The presented profiles, 


where the profile periods (peak-to-valley) correspond to a given feed value fz = 0.2 mm/tooth, testify to this. 


4 Gimadeev MR. Analiz zarubezhnogo opyta pri frezerovanii sferotsilindricheskim instrumentom. Khabarovsk: Informatsionnye tekhnologii 
XXI veka; 2019. P. 330-334. (In Russ.) 
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Fig. 10. Profilograms of surfaces after milling at different inclination angles y: 
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a— y= 10; direction f,; b — y = 10; direction a,; c — y = 40; direction f,; d—— y = 40; direction a, 


Photomicrographic images of the treated surface have an obvious hexagonal shape (Fig. 11). The quality of the 
machined surface in the direction of the side pitch a, may decrease due to the cutting tool retraction and an increase in 
the amplitude of the technological system vibrations under milling. This is validated by the results of vibroacoustic 
diagnostics. 125 
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b) 
Fig. 11. Milling scheme and photo of the treated surface at: a— y = 10; b— y =40 


We emphasize once again that the experimental results correlate with the data obtained by many authors, and are 
consistent with the theories of elasticity, shaping, material cutting, and mechanical engineering technology. 

Discussion and Conclusions. The influence of the inclination angle of the treated surface on the cutting forces, 
vibrations and microrelief parameters was studied. A model of cutting force and vibration was created, including 
kinematic and geometric parameters. The study has shown that quantitatively and qualitatively the cutting forces and 
vibrations of the cutting tool depend on the angle of inclination of the surface. 

An analysis of the experimental data made it possible to identify a trend towards a decrease in the amplitude 
parameters of roughness by 10—20 % with an increase in the angle of inclination to 30-40°. The quality of the treated 
surface was reduced due to the elastic retraction of the spherocylindrical tool and an increase in the amplitude of 
vibrations of the technological system. 

It was established that the forces grew in direction a,(X) and fz(Y) at y > 40°. This can be explained by the 
distribution of the cutting force components along the cutting edge, which depends on the inclination of the surface. 

The solutions to the problems of monitoring and analysis of roughness parameters obtained by the authors due to 
acoustic diagnostics can significantly reduce the volume of experimental studies. This approach makes it possible to 
significantly reduce the cost of material resources and time for mechanical processing of products. 

The experimental results are of practical interest and make the information about the formation of microrelief, 
including the roughness parameters under mechanical cutting, more accurate. A more detailed study of the proposed 


method of acoustic diagnostics is needed. 
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